SU(N) Evolution of a Frustrated Spin Ladder 
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Recent studies indicate that the weakly coupled, J± <C J, spin-1/2 Heisenberg antiferromagnet 
with next nearest neighbor frustration, J x , supports massive spinons for J x = J_l/2. The straight- 
forward SU(N) generalization of the low energy ladder Hamiltonian yields two independent SU(N) 
Thirring models with N — 1 multiplets of massive "spinon" excitations. We study the evolution of 
the complete set of low energy dynamical structure factors using form factors. Those corresponding 
to the smooth (staggered) magnetizations are qualitatively different (the same) in the N = 2 and 
N > 2 cases. The absence of single-particle peaks preserves the notion of spinons stabilized by 
frustration. In contrast to the ladder, we note that the N — > oo limit of the four chain model is not 
a trivial free theory. 

PACS numbers: 71.10.Pm, 75.10.Jm, 75.10.Pq 



I. INTRODUCTION 

Frustrated quantum antifcrromagncts are a source of 
considerable theoretical and experimental attention - 
see for example reference Ql Their characteristics include 
enhanced classical ground state degeneracies and the sup- 
pression of long-range Neel order. In addition to their in- 
trinsic interest, their prominance is fueled by the high-T c 
superconducting cuprates, where hole doping frustrates, 
and ultimately destroys the long-range Neel order of the 
parent compounds — see for example This motivates 
the quest for simple models of frustrated quantum mag- 
nets, and a detailed understanding of their properties. 

Important examples include nearest neighbor antifer- 
romagnets on frustrated lattices, such as the triangular, 3 
pyrochlorc, and Kagoma^ lattices, and further neigh- 
bor models on regular lattices. The second variety 
embraces frustrated chains^ and ladders^Si^ the pla- 
nar pyrochlorej2ii2iii and the square lattice antiferro- 
magnet with next-nearest neighbor interactions. In- 
deed, the latter model was suggested by Anderson in 
his influential worki^ on La2Cu04, as a means to re- 
alize his "resonating-valence-bond" (RVB) or "quantum 
spin-liquid" state. With isotropic nearest neighbor ex- 
change, J\, this is often referred to as the J\ — Ji model 
- for an introduction to spin-liquids see Chapter 6 of 
the book by Fradkin^ Other examples include multispin 
exchange models, and those of dimers^ Although enor- 
mous progress continues to be made, frustrated quantum 
magnetism remains theoretically challenging. In general 
one must resort to 1/S or 1/N expansions, numerical sim- 
ulations, or other approximation schemes — see for ex- 
ample reference 1 15i 

Building on the work of reference 0, Nersesyan and 
Tsvelik have made considerable advances in the so-called 
confederate flag models This is an anisotropic version 
of the much studied J\ — J% model, in which the nearest 
neighbor exchange has a strongly prefered chain direc- 
tion — see figure ^ The limit J x < J± <C J may be 
viewed as a collection of weakly coupled, but neverthe- 
less interacting chains, and field theory methods may be 
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FIG. 1: 2D Heisenberg antiferromagnet with next nearest 
neighbour frustration J x < J± <C J- The strongly relevant 
interchain interaction between staggered magnetizations van- 
ishes for J x — J±/2 and renders deconfined spinons. 



employed. In general, the massless spinons of the spin- 
1/2 chainii are confined by the interchain interactions. 
However, along the line, J x = J±/2, massive spinons 
emerge in pairs, as the elementary spin excitations of the 
coupled systemi2ii& In general they are neither bosons 
nor fermions, but have momentum dependent scattering. 
There have been many speculations about the existence 
of such excitations in two-dimensional frustrated antifer- 
romagnets, and their possible role in high-T c . The de- 
velopments of reference 1 16l deserve further investigation. 

In this paper we return to an SUfN) generalization of 
the ladder introduced in reference Our motivation is 
twofold: firstly, the large-N approach is known to miss 
qualitative features in this case^ and we wish to track 
its evolution in detail. Largc-N results will be important 
in two-dimensions, and we hope to gain expertise in all 
the solvable cases. Secondly, we calculate the dynamical 
structure factors of the staggered magnetizations. These 
involve correlation functions of interacting WZNW fields, 
and their evaluation beyond the ladder, is a highly chal- 
lenging and open problem^£*i£ 

The layout of this paper is as follows: in section HT1 we 
reacquaint the reader with the spin-1/2 model, and it's 
mapping on to two different "parity" sectorsi2ii& We in- 
troduce the SU(N) variant of the low energy action and 
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comment on this choice of generalization. We emphasize 
that this treatment is not the same as replacing lattice 
spins by SU(N) generators^ we expand on this in ap- 
pendix ^ where we comment on the connection between 
filling and SU(N) "spin" representations. In section ITTT1 
we calculate the dynamical susceptibilities corresponding 
to the uniform and staggered magnetizations. We con- 
clude in section llVl with results for the four chain model. 
In appendix IbI we discuss in detail, the excitations, scat- 
tering and form factors of the SU(N) Thirring model. 
We hope this may be of some assitance to the unfamiliar 
reader. 



II. MODEL 

In this section we reacquaint the reader with the spin- 
1/2 confederate flag model model, and it's mapping on 
to two different "parity" sectors^i^ we shall specialize to 
the ladder in due course. Consider a Hciscnbcrg antifer- 
romagnet on a two-dimensional square lattice (of spac- 
ing ao) with next nearest neighbor exchange interaction 
< J x < J± <S J as depicted in figure ^ 

L 

i— 1 n 

+ Jx (Si jTl • Sj_|_i jn+ i + S,- in _|_i • Sj_|_i jn )] . (1) 

It is well established, that the low energy dynamics of a 
single spin-1/2 (isotropic) Heisenberg chain 

Hi° = JSi,n ' Si,n+1, (2) 
n 

are described by the su(2)i Wess-Zumi no-Nov ikov- 
Witten (WZNW) model^2i for a review see l22l23i This 
WZNW model has conserved currents J = L\Jt a pCp 
and J = lZl,.t a plZp, which generate the su(2)i Kac- 
Moody current algebra, and the Hamiltonian density, 
H = J dxTl, may be written in the following (Sugawara) 
form: 

■H) U = Xhr (: J; .1, : (3) 

Here v is the spin velocity, Af is a normalization con- 
stant, and the ellipsis stand for less relevant operators. 
We replace the perturbing lattice spin operators by their 
continuous, slowly varying, uniform and staggered com- 
ponents: 

S^/ao -> S<(aO - M,(x) + (-l)"N;(x), (4) 

where x = nao measures the distance along chain i. Ne- 
glecting oscillatory and derivative terms, Hamiltonian 
becomes H = J dxH, where 

TV 
i=l 

+ 2J x a (Mi ■ Mi+i - N t ■ N. 1+1 ) . (5) 



In terms of the currents, Mj = 3{ + J;, the Hamiltonian 
density (JSJ may be written&i& 

N 
i=l 

+ A,\, • N , . ! • ••• (6) 

where 

Ai = (Jj_ + 2J x )a , \ 2 = (J±-2Jx)a Q . (7) 

In particular, for J x = Ji/2, the strongly relevant inter- 
chain coupling, A2, between the staggered magnetizations 
vanishesi£*i& Setting J x = J_l/2, and neglecting velocity 
renormalizing terms, the Hamiltonian splits into two in- 
dependent pieces, or "parity" sectors&±& 

H = H++H-, (8) 

where 

U+ = y^AfHv (J 2i ■ 3 2 i + Jat+i • Jgj+i) 

i 

+AiJ 2l ■ Jm+i, (9) 

and 7i_ is obtained from 7i + by the (parity) transfor- 
mation J J. In the sector of positive parity, the even 
(odd) chains carry left (right) moving fields; in the sector 
of negative parity the situation is reversed (see figure 
Equivalently, 7i+ and 7Y_ are interchanged under a shift 
by ao transverse to the chains. Specializing to the ladder: 

n+ = Nhv (Ji • Ji + j„ • j„) + Ax(Ji • j„), (10) 

where we label the chains by Roman numerals to avoid 
subsequent confusion with space-time indices. The 
Hamiltonian I|1U|) may be brought into a more familiar 
form by introducing a spinor, the left component of which 
resides on one chain and the right resides on the other: 

In terms of this spinor, the Hamiltonian becomes 

H+ =Afhv(J+- J++J+-J + ) + Ai(J + - J+) (12) 

and similarly for 7i_. (Equivalently one may perform 
the chiral interchange Ji «-> Jn on the original Hamilto- 
nians.) We see that H+ is nothing but an SU(2) Thirring 
model. That is to say, the frustrated ladder may be re- 
formulated as the sum of two decoupled SU(2) Thirring 
models, labelled by their parityi&i& We emphasize that 
each of these decoupled models captures the behavior of 
the coupled ladder, as highlighted in JTTJ, and not just a 
single chain. In particular, the elementary excitations of 
the ladder are those of the SU(2) Thirring model, namely 
massive spinons. These correspond to domain walls sep- 
arating regions of different spontaneous dimcrization. 8 
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FIG. 2: The ladder Hamiltonian is the sum of two indepen- 
dent SU(N) Thirring models: Ti = TL+ + H-. In the sector 
of positive "parity" the even (odd) chains carry left (right) 
moving fields. The sector of negative parity is obtained by 
reversing the arrows. Excitations of the ladder carry this ± 
index and may be produced in both sectors. 



In this paper we straightfowardly replace the SU(2) 
currents by SU(N) currents, as suggested in reference^. 
In each parity sector, the Hamiltonian becomes that of 
the SU(N) Thirring model with N — 1 multiplets of mas- 
sive spinons (see appendix lB|l . This is the simplest gen- 
eralization which retains spinon excitations and parity 
sectors. We note that the alternative strategy of replac- 
ing lattice spins by SU(N) generators leads to problems at 
the outset iiS As we discuss in appendix [S] the represen- 
tation of the generators translates into the filling of the 
corresponding electronic model. For the critical SU(N) 
Heisenberg model, with spins in the lowest fundamen- 
tal representation, the corresponding Hubbard model has 
one electron per site^ The corresponding "spin" den- 
sity has harmonics at multiples of 2k f = 27r/Nao due 
to all the fundamental primaries of the su(N)i WZNW 
model. In this case, the simple fine tuning condition, 
,/ x = Ji/2, does not remove all relevant perturbations. 19 
The absence of such terms is crucial for spinons in the 
confederate flag model, and such a generalization would 
be inappropriate. Attempts to reinstate the condition 
of half-filling with Hubbard chains or the alternating 
N eg) N magnet (q = N 2 quantum Potts model) also 
lead to difficulties; for N > 2 they are massive and 
dimerized2£i2§*2L2& and have little in common with the 
UV limit of decoupled spin-1/2 chains. Since our inter- 
est in these generalized models stems from the spinon 
physics of the confederate flag model, we confine our- 
selves to the simple minded extension of the low energy 
action. We study the SU(N) evolution of the original 
operators, and retain the terms smooth and staggered 
magnetizations for these fields. 

In the next section, we shall compute the dynamical 
structure factors of the generalized model. These are a 
direct probe of the elementary excitations. 



III. DYNAMICAL STRUCTURE FACTOR 

In this section we compute the dynamical structure fac- 
tor (as may be seen by neutrons) for momentum transfers 
close to the "soft modes" at and tt. This is nothing but 
a Fourier transform of the spin-spin correlation functions: 



S(u>, q, qj_) cx Imi 



dx 



( [Sf (t, x) ± S&(t, x), Sf (0, 0) ± S£(0, 0)] ). (13) 

The plus (minus) sign corresponds to q± = (q± = tt), 
and 6 ensures convergence of the temporal integral. The 
longitudinal momentum transfers in the vicinity of q = 
(q = tt) probe the smooth (staggered) components of the 
spin operators. The task is to relate the spin operators 
entering (jl3|) to the operators of the Thirring models, 
and to evaluate their matrix elements. 



A. Smooth Components 

The smooth component of the sum of the chain spin 
densities may be expressed in terms of the two Thirring 
models as follows: 



Si + Sn I smooth 



Ji + Ji + Ju + Jn 

Jo,+ + Jo,- 



where Jo,+ = Ji + Jn (Jo,- = Ji + Jii) is the temporal 
component of the Thirring current in the model of pos- 
itive (negative) parity. Simply put, the structure factor 
S(w,q ~ 0,0) of the frustrated ladder, may be obtained 
from the correlators of Jq in the SU(N) Thirring model. 



S(uj, q 



dx / 

-OC J 



£± i{uj+i8)t — ivqx 



[j >(t,z),Jo>(o,o)]; 



(14) 



where the summation is over parity sectors. The ele- 
mentary excitations of the SU(N) Thirring model consist 
of N — 1 multiplets of massive particles, corresponding to 
the fundamental representations of SU(N). The length of 
the Young tableau is termed the "rank" of the particle^ 9 - 
and their masses are given by (|B1|) . It is convenient to 
move to a basis of such particles and to parametrize their 
energy and momentum in terms of rapidity: 



E{ = rrij cosh 8j Pj = nii sinh ( 



(15) 



One may now insert a complete set of states between the 
current operators in Ijl4(l 



1^ 1^ I (9^2 n ) l y " 



?1=0 £; 



(27r) 2 n 



Me n ... ei ei -W0i.. 



(16) 

where the ti are the internal (or isotopic) indices carried 
by the members of each multiplct. Using 

<-^(0[...9' n \O(t,x)\e n ...9 1 ) en ... ei 

^■■■<{6[...e' n \o(o,o)\e n ...e 1 )z n ... ei (17) 

one obtains 

S(u, g~ 0,0) oc 
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S(vq — J2j m j sinh 9 j) 



S(vq + ^ Til j sinh 9j 



J2j m j cosh 6j + iS u> + J2j m j cosh 9j 



(18) 



where Fja(9 l . . . n ) ei ... 6n is a multiparticle form factor 
of the temporal Thirring current: 



Fjg(9 1 . . .9 n ) tl .. 



\J${0,0)\9 n ...9 



l/6„...ei 



(19) 



The dominant contributions to (|18fl come from the states 
with the lowest mass. In the case at hand these are two 
particle states of the (rank-1) fundamental □, and it's 
(rank-N-1) conjugate □. In particular, the current oper- 
ator couples to the adjoint representation occuring in the 
SU(N) tensor product □ ® □; for N = 2, □ is □. As we 
discuss in appendix iBl this form factor is 




2 / 2 

s /m 



FIG. 3: Exact dynamical structure factor S(u>, q ~ 0, 0) for 
N = 2 (solid), N = 3 (dashed) and N = oo (dotte d) with 
fixed q and arbitrary normalization. The y/ s 2 — Am 2 thresh- 
old behavior for the N = 2 physical case, is replaced by a 
1/Vs 2 — 4m 2 divergence for all N > 2. 



Fj {9i,9 2 ) a ,Q <*■ sinh 
where 



/ad? ( 



(20) 



/adjVi 



exp 




2 exp(a;/N) sinh(x/N) shr 2 (:ri/27r) 
x sinh 2 x 



(21) 



and 9 = in — 9; see equations <|B24I) and 1B27|) . We have 
supressed the isotopic and component information in 121|) 
and concentrated solely on the rapidity dependence. In- 
serting this into (|18fl and performing the 9 integrations 
one obtains 



5 , l/a D d?[2g(s)]| 2 (22) 



s 3 y / s 2 — Am' 2 



where s 2 



2 2 

v q , 



4rn 2 < s 2 < 



6*(s) = arcosh(s/2m) and 

16to 2 N = 2, 

9to 2 N = 3, 

16ra 2 cos 2 (7r/N) N > 3. 



(23) 



This result is plotted in figure [3] and is exact, provided 
(|23|1 is fulfilled. For larger energy transfers there are 
small corrections due to higher mass states; the upper 
thresholds correspond to four rank-1 solitons, three rank- 
1 (or rank-2) solitons, and a rank-2 bound state and its 
conjugate, respectively. In particular, there are no single- 
particle bound states appearing below the gap; the ele- 
mentary Thirring excitations correspond to fundamental 
SU(N) representations, and do not couple to the current 
directly, which spans the adjoint. 

The result (|22|) interpolates between two known limits. 
For N = 2, it coincides with equation (34) of reference H, 
and in the limit N = oo, where (|21() tends to unity, we re- 
cover the result for free massive fermionsA^ In particular, 
the (9 = 0) threshold behavior of (|21|l is quite instruc- 
tive: for N = 2 it vanishes like sinh(6>/2), as may be seen 



from (|B32|) , whereas it is finite and non- vanishing for any 
N > 2. As a result, the structure factor l|22|l vanishes as 
— 4m 2 in the physical case of N = 2, but diverges 
as l/Vs 2 — 4m 2 for any N > 2 — see figure 01 Solely on 
the basis of the N = 2 and N — oo limits)^ one might 
have expected the threshold to get steeper and narrower 
with increasing N, but to remain qualitatively correct for 
N < oo. The actual evolution, and the departure even for 
N = 3, is a sobering example of how SU(N) treatments 
may miss simple features over the entire range of N. 

Likewise, the smooth component of the difference of 
the chain spin densities may be expressed in terms of the 
two Thirring models as follows: 



*II smooth 



Ji + Ji — Jn 
Jl,+ — Jl,- 



where Ji, + = Ji — Jn (Ji,- = Ji Jn) is the spatial 
component of the Thirring current in the model of pos- 
itive (negative) parity. Simply put, the structure factor 
S(w, q ~ 0, 7r) of the frustrated ladder, may be obtained 
from the correlators of Ji in the SU(N) Thirring model. 
The corresponding form factor is given by (|B25|) : 



Fj 1 {9 1 ,9 2 ) d q oc to cosh 



/ a D d °(M. (24) 



We obtain 

S(w, q ~ 0, 7r) oc 



2 2 
77l Z UJ 



|/ a D d °[20( S )]| 2 . (25) 



s 3 \/ s 2 — 4m 2 



Once again, this result interpolates between the known 
N = 2 and N = oo results^ and the SU(N) approach 
leads to qualitatively incorrect results over the entire 
range of N > 2. 



B. Staggered Components 

We denote the staggered component of the spin on 
chain I, Sx(t, x)| s t a gg. by Ni(t, x). In the UV limit (cor- 
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responding to decoupled chains and m = 0) N(i, x) is a 
spinlcss su(N)i primary field with (full) scaling dimension 
A = f — f /N. For the ladder we propose the following for- 
mula for the long distance asymptotics of the real space 
correlation functions: 

flNi(t, x) ± N n (*, a;)].[Ni(0, 0) ± N„(0, 0)]) 

cx (Ni(t,x).Ni(0,0)) ± (Ni(f,aO.N n (0,0)) 
cxTO 2A [i4Ti(mr)±A^(mr)] + --- (26) 



where r = yz^ = Vx 2 ~ t 2 (v = 1) and K v (x) is Mac- 
donald's function^ also known as the modified Bcsscl 
function of the third kindi^i The dots stand for more 
rapidly decaying terms. In order to get a feel for this re- 
sult we begin by studying a few limits. In the limit N — > 
oo, A — > 1, each parity sector reduces to non- interacting 
massive fermions. More specifically, N, may be replaced 
by the fermion bilinear C\ a t a pHi t p + Ti\ a t a pCi./3 and 
one obtains 

(NlNj) cx {4d){TliRi) (27) 
(Ni.N„) oc (4Hn){Hi£h) (28) 
with the usual massive Dirac fermion correlators: 



= 2mK (mr) 



(29) 
(30) 



— see for example chapter 13 of the book [32. In equa- 
tions lj2"?|) and we see quite clearly that the corre- 
lators of staggered magnetizations are products of corre- 
lators from the sectors of different parityi&ifi by defini- 
tion the left and right moving fields on a given chain 
belong to different sectors. In coupling to the stag- 
gered magnetizations, the solitons are still created in 
pairs, but belong to different sectorsAi& In a given sec- 
tor (i.e. Thirring model) we thus require the matrix el- 
ements of single- soliton creation operators. The matrix 
elements of such operators have only recently become 
available iSiSiS The free fermions appearing in (|27|l and 
(|2"5|) for N — > oo, are replaced by chiral fields L S ,R S} 
which are non-local single soliton creation operators and 
carry the Lorentz spin, ±A/2, of a Thirring soliton$2£i2£ 
we take the plus (minus) sign for left (right) movers. 
These chiral fields are the components of an (interacting) 
su(N)i primary field, and the Lorentz spin is nothing but 
the UV conformal dimension. The single-solition form 
factors of such operators are governed (upto normaliza- 
tion) solely by their Lorentz transformation properties: 



{0\C s \9)^ m ^ 2 e Ae / 2 



A/2 -A0/2 



, (31) 



and their two-point functions are now readily computed: 



^— Tmch6-\-ixmsh6 (32) 



/z\ A / 2 
y-j 2A'a(Wzz) 



m 



(33) 



-Tinch6-}-ixmsh6 



m 2K (mV zz) 



(34) 
(35) 



where z = t — ix and r = it. The results for (TZ\TZ S ) and 
(1Z\C S ) follow by interchanging z and z. In particular, 
equation (|33|l first appeared in the study of weakly cou- 
pled one-dimensional Mott insulators^ Replacing the 
correlators in (|27|) and l|28|l with these more general ex- 
pressions, the result l|26[) follows immediately. 

Further, the Macdonald function has the asymptotic 
expansion given by equation 9.7.2 of reference l38t 



K\ (mr) 



2mr 



1 



M-l (/i-l)(/i-9) 
8mr 2!(8mr) 2 



+ 



(36) 



where /i = 4A 2 . The leading term in Q36[) is independent 
of A, and at separations r S> 1/to, the interchain and 
intrachain correlations (amusingly) coincide: 



(N o (t,ar).N 6 (0,0)) 



,2A-1 



-2mr 



(37) 



Coupling the chains together not only generates exponen- 
tially decaying interchain correlations, but also modifies 
the l/r 2A behavior within the chains. 

Substituting l|26|) into the definition (|13|l and effecting 
the Fourier transforms we obtain the following structure 
factors: 



2A 



S(ui, q ~ 71", 0) oc 



[s + Vs 2 - 4,m 2 ] + (2 m 



2A 



S(ui, q ~ 71", 7r) oc 



[s + Vi 



Am 2 



4m 2 



2A 



(2m) 



2A 



3VS 



4m 2 



(38) 



(39) 



where s 
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- (q—ir) 2 . In deriving these expressions the 



reader may find the integral representations (|32|l and 134|) 
more convenient. At threshold, S(u>, q ~ ir, 0) diverges as 
1/Vs 2 — 4m 2 for all N, and we plot this behavior in figure 
01 the large s behavior is s~ 2//N . Similarly, at threshold, 
S(u>, q ~ 7r, 7r) tends to a constant for all N. In contrast 
to the magnetization correlators, we obtain qualitatively 
similar results over the entire range of N. 



IV. CONCLUSIONS 

In this paper we have studied the SU(N) evolution of a 
frustrated spin ladder&i& The dynamical structure fac- 
tors corresponding to the smooth (staggered) magnetiza- 
tions are shown to be qualitatively different (the same) in 
the N = 2 and N > 2 cases. A robust feature which sur- 
vives however, is the absence of coherent single-particle 
excitations at low energies. This is in stark contrast 
to the unfrustrated ladder, and reinforces the notion of 
spinous stabilized by frustration^^ 
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SV,q~ 7T,0) 



2 / 2 

/m 



FIG. 4: Exact dynamical structure factor S(u>, q ~ 7r, 0) for 
N = 2 (solid), N = 3 (d ashed) a nd N = oo (dotted). The 
threshold behavior is l/\/s 2 — m 2 for all N, and we have nor- 
malized accordingly. The SU(N) approach leads to qualita- 
tively similar results over the entire range of N. 



In closing we note that the N = oo limit of the two 
chain model is a free theory, but this is not so in general. 
In particular, each parity sector of the four chain model 
may be viewed as two (non-chiral) su(N)i WZNW models 
coupled by currents: 

H+ = su(N) 1 +su(N) 1 +AJ+J + . 

With transverse periodic boundary conditions the cur- 
rents J+ = Ji + Jin, J+ = Jn + Jiv generate su(N)2 
Kac-Moody algebras, and it is convenient to write iiS^ifi 



Ha 



Z N + (su(N) 2 + AJ+J H 



The Zn parafermions^i describe gapless non-magnetic 
excitations with central charge c = 2(N — 1)/(N + 2). 
They are unaffected by the interaction (due to the bound- 
ary conditions) and for N — * oo they reduce to two Gaus- 
sian models. The su(N)2 model on the other hand is 
rendered massive by the interaction, and is in fact inte- 
grable. The mass spectrum coincides with that of the 
(two chain) SU(N) Thirring model i|Bl|) . but the scatter- 
ing is notably different. The S-matrix can be extracted 
from a straightforward generalization of the Thermody- 
namic Bethe Ansatz (TBA) equations derived in refer- 
ence 0. At low temperatures (T <C Mi) the free energy 
of the perturbed su(N)k WZNW model is given by 



where * denotes convolution, the kernels C nm (X) and 
Aij(X) are given in reference and A = N0/27T. We 
extract the Bethe equations E = J22=i m ch9 a and 



exp(imL sh8 a ) 
Y[S o (0 a - 

b^a 



a 



where 



S (8) =exp{- f 



(1 + e -2|wk/N)2 



sh7r(l - l)w 
shnu) 



and 



e n (x) 



x — inn/N 
x + wrn/N 



£{x) = 



a Nz/2 



-Nx/2 



The rapdities X a and are distributed according to the 
A N_1 heirarchy, the details of which do not concern us 
here. Similar equations occur for the SU(N) invariant 
Thirring model (k = 1) but without the fi rapidities and 
with a different Sq(9). In the limit N — > oo one obtains 



exp(imLsh6 a ) = Y[S%°{9 ab ) 



b^a 

where S^°(8) = cxp (— i-| sign 6*) . This is to be con- 
trasted with the SU(N) invariant Thirring model where 
Sq 2 (6) = — 1. The absence of a simple N — > oo limit 
will be crucial for multiparticle form factors, and renders 
excitations with non-trivial statistics. In future publi- 
cations we hope to study these pertinent issues in more 
detail. Recent progress on spinon propagation in the four 
chain model may be found in the work of Smirnov and 
Tsveliki^ 
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j'=i J 
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ch6> ln[l + e" 



»(e)/ri 



where in this case k = 2. The excitation energies e 
(j = 1, . . . N - 1, n = 1, 2 . . .) satisfy 



Tln(l + e e " ) ( A )/ T ) 



ln(l 



5 n>h Mi ch(27rA/A0 



APPENDIX A: SPIN OPERATORS 

In this appendix we comment on the connection be- 
tween SU(N) spin representations and filling. At each lat- 
tice site (labelled by n) one may introduce the fermionic 
spin operators 



N 

£ 



L n,a''a/3 L n,P 



(Al) 
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where c and obey the canonical fermionic anticommu- 
tation relations 

{<4,a> c ™,/3} = 6 n , m 8 at p {c, c} = {c f , c 1 } = (A2) 

and the generators t a span the algebra su(N): [t a ,t & ] = 
if ab c t c . It is readily verified that spins on different sites 
commute, whereas those on the same site satisfy the 



su(N) algebra: [S%,S° 



f ab c S^. In the fun- 



damental representation, the generators are chosen to 
satisfy tr {t a t b ) = C5 ab , t a t a = C 2 I, with C = 1/2 and 
C 2 = (N 2 - 1)/2N; see appendix A.3 ofE3 

One may specify the su(N) representation on which 
spin operators S„ act by the relevant Young tableau — 
see for example l43t In particular, this fixes the value of 
the quadratic Casimir S 2 , and thus by equation JAIL 
constrains fermion occupation numbers. As we shall 
demonstrate, the constraint 



APPENDIX B: SU(N) THIRRING MODEL 

In this appendix we discuss the excitations, scattering 
matrices and form factors of the SU(N) Thirring (chi- 
ral Gross-Neveu) model. More details may be found in 
appendix A of Smirno^ and the literaturei 36 i 37 ' 46 i 47 ' 48 



1. Excitations 

The excitations of the SU(N) invariant Thirring (chiral 
Gross-Neveu) model are N — 1 multiplets of fundamental 
particles, corresponding to the N — 1 fundamental repre- 
sentations of SU(N). Their masses are given by 

M a = m S ^§; a = l,2,...,N-l, (Bl) 
sm7r/N 



^ 4i,a C n,a = h] Vn. 



(A3) 



and following Smirnov, we shall refer to the label a as 
the "rank" of the particle. 



corresponds to the vertical (i.e. antisymmetric) Young 
tableau of height h, as depicted in figure |S] The con- 



2. S-Matrices 



FIG. 5: Antisymmetric su(N) Young tableau of height h. 

straint (|A3(I fixes h electrons per site, and the permiss- 
ablc states to be of the form: 



,t 



JO), 



(A4) 



where € [1,---,N]. By virtue of the fermion anti- 
commutation relations 1A2L this may be viewed as a 
tensor of rank h, antisymmetric under the interchange 
of any pair of labels a; by the standard conventions for 
Young tableau 4 ^ this corresponds to a vertical diagram of 
h boxes. Moreover, it also follows from the anticommu- 
tation relations IA2L that there are N(N— 1) • • • (N — h + 
l)/h\ independent states of the form (|A4|I : this coincides 
with the dimension of the representation corresponding 
to the Young tableau of figEJ see §8.4 of Further, 
squaring equation (|A1|) and enforcing the constraint (|A3|) 
one obtains^ 



h(N 2 - h) h(l - h) 
2N + 2 ' 



(A5) 



This coincides with the quadratic Casimir of the su(N) 
Young tableau depicted in fig|SI 4 i see eq. 2.19 of0. e.g. 
for the fundamental □ of su(2) (h = 1,N = 2) one obtains 
S„ = 3/4, as appropriate for spin-1/2. 

Thus, equation (|A1|) supplemented by the constraint 
(IA3|I leads to spin operators S„ described by the Young 
tableau of figEJ 



The S-matrix describing the scattering of two (rank- 
1) fundamental particles in the SU(N) invariant Thirring 
(chiral Gross-Neveu) model is given bji 2 ^ 7 . 



S aa (9)\t 



(B2) 



where 8 = 6\ — 62, e € [1, ■ ■ • , N], and the S-matrix oper- 
ator acts on the two body Hilbcrt space □ ® □: 



(0)=S O (6) 



2m 
N 



(B3) 



X and 7^12 are the identity and permutation operator re- 
spectively with matrix elements 



e' 2 e[ (02#1 |^|#l#2) ei e 2 
e[ (020l|7 , 12|0102)eie 2 



r) 2 r) 1 



and 



S (6) 



2ni ' 



r(i 



N 



2ni 



(B4) 
(B5) 

(B6) 



See equation (11a) of |33 or appendix A (p. 182) of l29l 
e.g. for N = 2 this reduces to equation (6) of [2^. Using 
the decompositions 1 = + V { '\ and V12 = - 
V^~\ one may also write (JB31 in the form: 



s nn (6) = J2 s ? n wp (r) 

r 

= S (6) ( ?W + - n 



2tH > 
_N_ V (-) 

N j 



(B7) 
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where and "P^ - -* act on the symmetric and anti- 

symmetric representations occuring in the tensor prod- 
uct □ ® □; e.g. 3(8)3 = 6 + 3 in SU(3). Bound states 
correspond to poles of the S-matrix, with masses 



m b = 



m\ + 2m\m,i cosh(#i 2 ). 



(B8) 



Since T(z) is free of zeros, and exhibits simple poles at 
2 = 0, —1, —2, . . .j2i it follows that (|B7|) has a single sim- 
ple pole at 9 = 2iri /N occuring within the physical strip, 
< 9 < ni. This yields the bound state mass of the sec- 
ond fundamental particle, M 2 = m sin(27r/N)/ sin(7r/N), 
as given by equation I|B1|) . 

The S-matrix describing the scattering of a (rank-1) 
fundamental particle off its conjugate (rank-N — 1) may 
be obtained from l|B7|l by the crossing transformation: 



(B9) 



S D °(9) = C n S aa (iir-9)C a 



where is the conjugation operator on □. Utilizing 
■P (0) = C a V12 C n /N, and X = "P (adj) + V (a) one obtains 



= -Si(6) ( V 



)( ad j) 



>(0) 



where 



Si(0) 



rri + J-) rfi - i - -M 



rfi - J-) rfi 



' 2-iri I 



(BIO) 



(Bll) 



and p( ad J) and P (0) act on the adjoint and singlet rep- 
resentations occuring in the tensor product □ ® □; e.g. 
3(8)3 = 8 + 1 in SU(3). In particular, for N = 2, (JbTJ) and 
(|B10J| coincide (up to sign) as expected from the identi- 
fication of □ and □ in SU(2). Moreover, equation l|B7l) 
is also in agreement with equation (1.6a) of reference 149. 
The S-matrix l|B10[) has a pole at 9 = ni — 2iri /N occur- 
ing within the physical strip, < 9 < iri. This is a cross 
channel pole. 

For the purpose of calculating form factors in section 
IB 31 it proves useful to have the S-matrices in an inte- 
gral form. Taking the logarithm of (|B6|) and I|B11(I and 
employing 



lnr(z) 



one obtains 



dt 

T 



(z-l)c 



(B12) 



S Q (0)=ex P j^ dxf a (x) sinhf^ J (B13) 



where 



= 2exp(,/N)sinh[,(l-l/N )]; 



hi*) 



x sinh x 
2exp(a;/N) sinh(a;/N) 
a; sinh x 



(B15) 



In particular, the su(2) Thirring S-matrix coincides with 
the sine-Gordon S-matrix with 1 = 8ir& 



S a (9) — ? exp ^ % 



exp(-7TK/2) . 

uk — sin k9 

k cosli7Tfi;/2 



(B16) 



Form Factors 



In the previous paragraphs, we have discussed the ele- 
mentary excitations of the SU(N) Thirring model. They 
are massive particles labeled by their rapdities, 9i, and 
carrying quantum numbers or isotopic indices, e^. In or- 
der to compute correlation functions and dynamical sus- 
ceptibilities, we will need the matrix elements of various 
physical operators, 0, between the vacuum and the (low- 
est) multiparticle excited states. Such matrix elements 
are termed form factors, and their computation is an im- 
portant enterprise; see for example l29l As is discussed in 
Ch. 1 of Smirnov's book^ the two-particle form factors 



F O (6l,0 2 )e u e 2 = (O|0(O,O)|0 2 0i) £2 , ei 



(B17) 



satisfy a matrix (Ricmann-Hilbert) problem, also known 
as Watson's equations: 

F(6 U 9 2 + 27u) ei , £2 = F{9 1 ,92)^Stri{9 12 ). (B18) 

This equation may be diagonalized to yield the simpler 
scalar problcm(s) 



F(9 1 ,9 2 + 2m) =F(6 u e 2 )S(e 12 ), 



(B19) 



where S(9) are the S-matrix eigenvalues. In particular, 
the Thirring current operator (with N 2 — 1 compo- 
nents) couples to the adjoint representation occuring in 
the tensor product □ <8> □; the relevant eigenvalue is 

S(0)=S°J(0) = -Si(0). (B20) 

Another constraint on the form factors comes form 
Lorentz invariance. Under a Lorentz boost, correspond- 
ing to a simultaneous shift of all rapidities by A, the 
two-particle form factor of an operator O of spin s satis- 
fies: 



F {6t + A, 9 2 + A) = e sA F (9 1 ,6 2 ) 



(B21) 



In particular, the left (right) component of the Thirring 
current has spin ,s = +1 (s = —1) and one obtains: 

FftiOuO*) (x e + ^ +e ^ 2 f°f(9 12 ), (B22) 
F j% (9 1 ,9 2 ) « e - (ei+e2)/2 / a D d°(M. (B23) 



Note that f^(9i 2 ) is a function of 9\ — 9 2 , and is thus 
Lorentz invariant. The form factors corresponding to the 
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temporal and spatial components of the current may be 
written: 



/ 9 

Fjg (6*i , 62 ) <x m sinh I — 



1" u 2 \ 
2 I J&dj 



f^(9 12 ), (B24) 



F j? (6 1 ,6 2 ) <x mcosh ( -1±J. ) / a D d °(0 12 ). (B25) 



Substituting (|B20|) and either of ljB24|) and (|B25)l into 

'adj 

cUU/n o -\ -?□□/ 



(|B19|) . one obtains a constraint on 



/ a u d y(0-2«) = / a L d y(0)5 1 (^). 



(B26) 



Following the general arguments of Karowski and Weisa^ 
(equations 2.18 and 2.19) equation i|B26() may be solved 
by 



/°°(0)=cxp. 



[°° sin 2 (xg/27r) 

/ dxh(x) — 

/ sinh x 



(B27) 



where 6 = in — 9ml Expanding the denominator factors 
in powers of e~ 2x , and employing the identity 



£(j+!+ii) 



r(| 



+ Z- -M 



1 V 2 N 1-Ki) 



(B30) 



As may be seen most clearly from l|B29(l and l|B30|) , this 
form factor is free of poles in the physical strip < 9 < 
in, and Watson's minimal equations (in Karowski- Wcisz 
form) are explicitly satisfied: It is indeed, a minimal form 
factor. Expressions l|B27|l . I|B29|) and l|B30|) conform to 
the Karowski- Weisz normalization F(in) = 1, and have 
the asymptotic behaviour 



lim 



exp(±0/2N). 



(B31) 



In the limit N = 2, one may write l|B27|l in the form 



exp 



^2e-^sinh 7 x = te, 
x p — 7 



(B28) 



one obtains the equivalent representation: 



= n 



l,m=0 



1 + 1 + m 



1 ~ W + 1 + m 



\-± + l + m +J- 

h + i + m + ^- t 



2 N ' ' 27ri 



(B29) 



Application of Euler's Formula yields: 



/ a D d °w = n 



r(i-4 + 



r(i + o 



too 




ismh(9/2) x 

sin 2 (x0/27r) 
x sinh x 



[tanh(x/2) - 1] HB32) 



and expressions (|B24|) and (|B25|) coincide with the known 
results for the SU(2) invariant Thirring (or sine-Gordon) 
model; see equation (33) of Allen et aim or let £ — * 00 in 
the formula for / M (/3i, ffa) given on page 46 of SmirnoAi 2 ^ 
and note the different definition of the physical strip. 
In the limit N — * 00, the SU(N) Thirring model maps 
onto a theory of free massive fermions, as reflected in the 
explicit S-matrices. In this limit f^(0) -> 1, and (|B24|) 
and i|B25|l coincide with the free fermion form factors 
given in equation (108) of Smirnov^ 
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